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Abstract 

Recently, it has been demonstrated that an arbitrary hnear quantum stochas- 
tic system can be reahzed as a cascade connection of simple one degree of freedom 
quantum harmonic oscillators together with a direct interaction Hamiltonian which is 
bilinear in the canonical operators of the oscillators. However, from an experimental 
point of view, realizations by pure cascading, without a direct interaction Hamilto- 
nian, would be much simpler to implement and this raises the natural question of 
what class of linear quantum stochastic systems are realizable by cascading alone. 
This paper gives a precise characterization of this class of linear quantum stochastic 
systems and then it is proved that, in the weaker sense of transfer function realizabil- 
ity, all passive linear quantum stochastic systems belong to this class. A constructive 
example is given to show the transfer function realization of a two degrees of freedom 
passive linear quantum stochastic system by pure cascading. 

Keywords: Linear quantum stochastic systems, quantum system realization, quantum 
networks, quantum control, linear quantum optics 

1 Background and Motivation 

Recently, there has been interest in the literature on control of a linear quantum stochastic 
system with a controller which is a quantum system of the same type ^ El El S] , often 
referred to as "coherent-feedback control". The potential applications for linear quantum 
stochastic systems include quantum information processing and photonic signal processing. 
For instance, they can act as the coherent photonic circuitry subsystem in a cavity QED 
system, the latter system being realized by placing suitable atoms inside the optical cavities 
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in a linear quantum stochastic system. Cavity QED networks are of interest for quantum 
information processing (see, e.g., [5]), such as in the quantum internet [6], whilst the 
controller realized in [1] is an early sample application of linear quantum stochastic systems 
to photonic signal processing. 

The studies on coherent-feedback control naturally led to the consideration of the net- 
work synthesis problem for linear quantum stochastic systems [7], which may be viewed 
as a quantum analogue of the network synthesis problem for linear electrical systems [8]. 
Nurdin, James and Doherty ^ have shown that any linear quantum stochastic system 
can, in principle, be synthesized by a cascade of simple one degree of freedom harmonic 
oscillators together with a direct interaction Hamiltonian between the canonical operators 
of these oscillators. Alternative schemes have subsequently been proposed in [21 [10], but we 
note that [10] considers a weaker type of realizability than in [ZIIH], i-e., transfer function 
reahzability (cf. section |3]), and the results therein limited to a certain sub-class of linear 
quantum stochastic systems. 

From an experimental perspective, direct bilinear interaction Hamiltonians between 
independent harmonic oscillators are challenging to implement for systems that have more 
than just a few degrees of freedom and therefore it becomes important to investigate 
what kind of systems can be realized by a pure cascade connection. A key result of this 
paper is a necessary and sufficient condition for a linear quantum stochastic system to 
be realizable by only a cascade connection of one degree of freedom oscillators, without 
any direct interaction Hamiltonian. Moreover, we also show that the associated transfer 
functions of all passive linear quantum stochastic systems can always be realized by a 
cascade connection, proving in general the partial results of [lO] without the additional 
assumptions made therein. 

The organization of this paper is as follows. Section [2] sets up the notations and gives a 
brief overview of linear quantum stochastic systems. Section |3] defines the synthesis problem 
and discusses the notions of strict realizability and transfer function realizability. Section |4] 
derives a necessary and sufficient condition for a linear quantum system to be realizable by 
a pure cascade connection of one degree of freedom quantum harmonic oscillators, in both 
the strict and transfer function sense of realizability. Section [5] then introduces the class 
of passive linear quantum systems and proves that all such systems are transfer functions 
realizable by a pure cascade connection. Finally, section |6] offers some conclusions of this 
paper. 

2 Preliminaries 
2.1 Notation 

We shall use the following notations: i = * denotes the adjoint of a linear operator as 

well as the conjugate of a complex number, li A = [ajk] then A* = [a*^], and = (A*)^, 
where denotes matrix transposition. 3?{y4} = {A + A*)/2 and Q'{v4} = |r(A — A*), 
and denote the identity matrix by / whenever its size can be inferred from context and 
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use In to denote an n x n identity matrix. Similarly, denotes a matrix with zero entries 
whose dimensions can be determined from context. diag(Mi, M2, . . . , M„) denotes a block 
diagonal matrix with square matrices Mi, M2, . . . , M„ on its diagonal block, and diag„(M) 
a block diagonal matrix with the square matrix M appearing on its diagonal blocks n 
times. 

2.2 The class of linear quantum stochastic systems 

In this paper, we will be concerned with a class of quantum stochastic models of open 
(i.e., quantum systems that can interact with an environment) Markov quantum systems 
that are widely used and are standard in quantum optics. Such models have been in the 
physics and mathematical physics literature since the 1980's, see, e.g., [HI [121 HSl El US]- 
In particular, we focus on the special sub-class of linear quantum stochastic models, see, 
e.g., [151 section 6.6], [HI sections 3, 3.4.3, 5.3, chapters 7 and 10], [16, section 4], [T71 
section 5] , [HI El [3l [3, UHl HI EDI El] . These linear quantum stochastic models describe such 
quantum optical devices as optical cavities [221 section 5.3.6] [231 chapter 7], linear quantum 
amplifiers [HI chapter 7], and finite bandwidth squeezers [HI chapter 10]. Following the 
terminology in [21 |3l [7] , we shall refer to this class of models as linear quantum stochastic 
systems. 

Suppose we have n independent quantum harmonic oscillators labelled 1, . . . ,n. Each 
oscillator j has position and momentum operators qj and Pj, respectively. The position 
and momentum operators satisfy the canonical commutation relations [qj,Pk] = 2i5jfc, 
fej^fc] = 0, and \pj,Pk] = 0, where 6jk denotes the Kronecker delta that takes on the 
value 1 only if j = k, but is otherwise 0. Equivalently, we may describe them in terms of 
the 2n annihilation and creation operators ai, a^, 02, a^, ■ ■ ■ , an, a*, with aj = {qj + ipj)/2, 
satisfying the canonical commutation relations [aj,^^] = Sjk, [aj,a/c] = and [0^,0^] = 0. 
The independent oscillators can be coupled to one or more external independent quantum 
fields, say m of them. In a Markov quantum system, the m independent fields are essentially 
quantum noises modelled by bosonic annihilation field operators Ai(t) , A2(t) , . . . , Am(t) 
that can be defined on a separate Fock space (over L^(]R)) for each field operator [HI [HI El]. 
For each Aj(t) there is a corresponding creation field operator A*(t) that is defined on the 
same Fock space and is the operator adjoint of Aj{t), i.e., A*it) = Aj{t)*. The field 
operators are adapted quantum stochastic processes with forward differentials dAj{t) = 
Aj{t + dt) - Aj{t) and dA*{t) = A*{t + dt) - A*{t) that have the quantum Ito products 
[HI [HI El: 

dAj{t)dAk{ty = 6jkdt; dA*{t)dAk{t) = 0; dAj{t)dAk{t) = 0; 
dA*{t)dAl{t) = 0; dAk{t)dt = 0; dAl{t)dt = 0. 

More informally, as in the quantum Langevin formalism, we can express Aj{t) = j^rjj{s)ds 
and A*{t) = Jp r]*{s)ds, where r]j{t) for j = 1, . . . , m are independent quantum white noise 
processes satisfying the informal commutation relations [^7j(s), ?7^(t)] = SjkS{t — s) and 
[?7j(s), ?7A;(i)] = [Vji^y :Vki't)] = 0) where ri*(t) = rjjit)*, and 6(t) denotes the Dirac delta 
function. 
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Let us collect the position and momentum operators in the column vector x defined 
as X = {qi,Pi,q2,P2, ■ ■ ■ ,fln,Pn)'^ ■ Note that in terms of x, we may write the canoni- 
cal commutation relations as xx'^ — {xx'^Y' = 2i0 with = diag„(J). We take the 
composite system of n quantum harmonic oscillators to have a quadratic Hamiltonian 
H given by H = ^x'^Rx, where i? is a real 2n x 2n symmetric matrix. The oscilla- 
tors are coupled to the quantum field m via the informal singular interaction Hamil- 
tonian Hm = i(L^r7;;(t) - L*^Vm{t)) [El H, where L„ = K^x with G C^^^n jg 
a linear coupling operator of the oscillator position and momentum operators to r]m{t). 
Collect the coupling operators Li, L2, . . . , together in one linear coupling vector L = 
{Li,L2, . . . ,Lm)^ = Kx, with K = [ ■ ■ ■ K'^ ]'^, and the field operators to- 

gether as A{t) = {Ai(t) , A2{t) , . . . ,Am{t))'^. Then the joint evolution of the oscillators 
and the quantum fields is given by a unitary adapted process U (t) satisfying the Hudson- 
Parthasarathy quantum stochastic differential equation (QSDE) [TT ] IT ^ 1^: 

dU{t) = (tr((S - lfdA{t)) + dA{t)^L - L^SdAit) - {\H + ^L^ Ldt))U{t), 

where S G C™"^™ is a complex unitary matrix (i.e., S'^S = SS'' = I) called the scattering 
matrix, and A(t) = [Ajk(t)]j^k=i,...,m- The processes Ajkit) for j, k = 1, . . . ,m are adapted 
quantum stochastic processes that are referred to as gauge processes, and the forward 
differentials dAjk(t) = Aj^it + dt) — Ajkit) j, k = 1, . . . ,m have the quantum Ito products: 

dAjk it)dAj,k' (t) = 5kj'dAjk' (t) , dAj {t)dAki (t) = 5jkdAi (t) , dAjkdAi {t)* = 5kidA* (t) , 

with all other remaining cross products between dAj^it) and either of dt, dAfit) or dAl,{t) 
being zero. Informally, we may express Ajkit) = ri*{s)rik{s)ds. 

For any adapted processes V{t) and W(t) satisfying a quantum Ito stochastic differen- 
tial equation, we have the quantum Ito rule d{V{t)W{t)) = V{t)dW{t) + {dV{t))W{t) + 
dV(t)dW(t). Using the quantum Ito rule and the quantum Ito products given above, as 
well as exploiting the canonical commutation relations between the operators in x, the 
Heisenberg evolution X{t) = U{t)*xU{t) of the canonical operators in the vector x satisfies 
the quantum stochastic differential equation, see [16], section 4], [T71 section 5], [2l[7j: 



dX{t) = d{U{t)*xU{t)) = AX{t)dt + B 



X(0) = X, 



dA{t) 
dA{t)* 

dY{t) = d{U{t)*A{t)U{t)) = Cx{t)dt + DdAit), (1) 



with A = 2e{R + ^{K^K}), B = 2i0[ -K^S K^S* ], C = K, and D = S, where 
Y{t) = (Fi(t), . . .,Ym{t))^ = U{t)*A{t)U{t) is a vector of output fields that resuhs from 
the interaction of the quantum harmonic oscillators and the incoming quantum fields A{t). 
Note that the dynamics of X{t) is linear, while Y{t) depends linearly on X{t) and A{t). 
We refer to n as the degrees of freedom of the oscillators. If = 1, we shall often refer to 
the linear quantum stochastic system as a one degree of freedom (open quantum harmonic) 
oscillator. 

Following [25j, we denote a linear quantum stochastic system with Hamiltonian H, 
coupling vector L and scattering matrix S simply as G = {S, L, H) or G = {S, Kx, \x'^Rx). 
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We also recall the concatenation product ffl and series product < for open Markov quantum 
systems [25] defined by Gi ffl G2 = idiag{Si, S2), (Lj , L^f , Hi + H2), and G2 < Gi = 
{S2S1, L2 + S2L1, Hi + H2 + '^{lIS2Li}). Since both products are associative, the products 
GiS G2S . . . S On and Gn < Gn-i < . . . <Gi are unambiguously defined. 

3 Synthesis of linear quantum stochastic systems 

The network synthesis problem for linear quantum stochastic systems can be stated (in a 
strict sense, as explained below) as the problem of how to systematically realize a given 
linear quantum stochastic system with a given fixed set of matrix parameters S, K, R from 
a bin of certain basic quantum optical components; see for details of these basic compo- 
nents. A particular solution was proposed to the synthesis problem, see [U Theorem 5.1]: 
Any linear quantum stochastic system with n degrees of freedom can be synthesized via a 
quantum network consisting of a cascade connection of n one degree of freedom harmonic 
oscillators together with a direct interaction Hamiltonian that is bilinear in the canonical 
operators of the oscillators. Partition R as R = [Rjk]j,k=i,...,n with Rjk G M^^^ and K as 
K = [ Ki K2 ■ ■ ■ Kn ] with Kk G C™^^. Then according to [7T, Theorem 5.1] a system 
G = {S, Kx, Rx) can be decomposed as G = {Gn<Gn-i<- ■ •<Gi)ffl(0, 0, if"'), where the 
GiS are (simpler) one degree of freedom open harmonic oscillators Gi = {Si, KiXi, ^xjRaXi) 
{xi = {qi,PiY) with parameter values specified by the theorem, and H'^ is a direct bilin- 
ear interaction Hamiltonian of the form H'^ = X]J=i Ylk=j+i^J {^jk ~ ^{-^I-^il^j ^k- 
The work [7j then shows how each of the Gj's can be synthesized from the bin of given 
components and how H"^ can be realized. However, in current practical experiments, im- 
plementation of H"^ can be challenging for systems that have more than just a few degrees 
of freedom. Therefore, it is of interest to characterize the class of systems that can be 
synthesized by pure cascade connection alone, that is, with H'^ = 0. 

As alluded to at the beginning of this section, we emphasize that considers a strict 
type of realization problem, that is, it deals with how to synthesize a given and fixed triplet 
{S, L = Kx, H = ^x'^Rx} that describes a linear quantum stochastic system G. This type 
of strict realizability is relevant, for instance, in cases where the internal dynamics X{t) 
may represent some (continuous time) quantum information processing algorithm and thus 
needs to be realized as given. However, for some linear quantum control problems such 
as robust disturbance attenuation j2j and LQG synthesis [3], the internal dynamics are 
inconsequential. In this case there is freedom to modify /transform these dynamics and 
what is important is the associated (classical) complex transfer function associated with 
the system matrices {A, B, C, -D)Q As is well known, a transfer function is invariant under a 

^As in [2], here we shall not define the transfer function of quantum systems, but associate to a 
quantum system G with system matrices {A, B, C, D) a classical, doubled-up [26], transfer function G{s) = 
[ — A)~^^B + dia.g{D,D'^). However, we also remark that G(s) can actually be properly 

interpreted as a genuine transfer function for the quantum system following [^Hl [HI [571 H] ) this being a 
common practice in the physics community via Fourier transform methods |23[ I14| . In any case, we are 
dealing with the same object G{s) and thus the particular interpretation attached to it becomes immaterial 
for our purpose. 
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similarity transformation of the system matrices {A,B,C,D) (yAV~^,VB,CV~^,D) 
for any invertible matrix V. However, for linear quantum systems the transformation 
matrix V for a similarity transformation is restricted in that it has to be a symplectic matrix: 
V is real and satisfies the condition VQV^ — ©. This ensures that the transformed variable 
Z{t) = VX{t) satisfies the required canonical commutation relations (CCR) of quantum 
mechanics: Z{t)Z{t)^ - {Z{t)Z{t)^)^ = 219, so the system remains physical. Note that 
the set of all symplectic matrices of a fixed dimension form a group and in particular 
V^^ is again a symplectic matrix. Such a similarity transformation in quaiituui systems 
corresponds to replacing G = {S,Kx,^x^Rx) with G' — {S,KV~^x,^x-'^V~^RV~^x). 
This motivates us to introduce the following definition: 

Definition 1 Let G — {S, Kx, ^x'^Rx) and G' — {S', K'x, \x'^R'x) he two linear quantum 
stochastic systems. Then G' is said to he transfer function equivalent to G or is a transfer 
function realization of G if S' ~ S and there exists a symplectic matrix V such that 
R' = V-^RV-\ K' = KV-^ (or, equivalently, R = V^R'V and K = K'V). G is then 
said to he transfer function realizable hy G' , and vice-versa. 

Renicirk 2 It is important to note that two transfer function equivalent systems G and 
G' will not necessarily generate the same input-output dynamics {A(t),Y(t)) for all t>0. 
This is hecause although they can have different parameters, they always have the same 
initial value X{0) = x, whilst for quantum systems x clearly cannot he zero due to the 
CCR condition. If the A matrix of G is Hurwitz then the input-output dynamics of G and 
G' converge in the limit t ^ oo. However, as remarked earlier, for some linear quantum 
control design ohjectives internal dynamics and initial conditions do not play an essential 
role, only the transfer function does. 



4 Conditions for realizability by a pure cascade con- 
nection 



In this section we state and prove a theorem that characterizes the class of linear quantum 
stochastic systems that can be realized simply by a cascade connection of one degree of 
freedom (open quantum harmonic) oscillators. Let us first introduce the following notation: 
Sk^i ^ Sk--- Sj+iSi for aU j < k, Sk^k = Sk and Sk^k+i = Im, and let Xi = {qi.Pif for 

„T\T 



'k'>^3 — ^fe ■ ■ ■ '^j+l'-'j 

i = 1, . . . , n SO that x 



(xf , . . . , x^) , where xx — {xx ) — 210. Moreover, we introduce 
the following terminology: A square matrix F is said to be lower 2x2 block triangular if 
it has a lower block triangular form when partitioned into 2x2 blocks: 



Fnl 



02x2 
-^22 



02x2 
02x2 



02x2 
02x2 



n2 



where Fjk, j < /c, is of dimension 2x2. We start with the following lemma: 
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Lemma 3 The cascade connection Gn<Gn-i <■ ■ ■ <iGi of one degree of freedom harmonic 
oscillators Gi = {Si, KiXi, ^xj RiXi) [i = l,...,n) realizes a linear quantum stochastic 
system G = {S,Kx, \x^ Rx) with S = Sn^i, K = [ Sn<^2Ki Sn<^zK2 ■ ■ ■ Kn ] , R = 
[Rij]i,j=i,...,n, where Rjj = Rj, Rkj = Q{KlSk<^j-^-iKj} whenever k < j and Rjk = R^j 
whenever j > k. In particular, R + '^{K'^K} is lower 2x2 block triangular. 



Proof. The proof proceeds along the hnes of the proof of |7l Theorem 5.1]. By the series 



product formula (cf. section 2.2) for the cascade of two one degree of freedom oscillators 

Gi = {Si,KiXi,^xjRiXi) and G2 
G2<G, 



(X^ ,x 



T\T 



{S2, K2X2, 1x2 R2X2)^ we get the oscillator G(2) 
(5'25'i, S2K1X1 + K2X2, |a;f + ^X2R2X2 + xlQ{KlS2Ki}xi). Letting X(2) 



1 5 -^2 . 



(2) 



s. 



the latter may be compactly written as: G(2) = (5'(2), K(^2)X{2), | ^(2) -^(2) 2^(2)) with 

Rl '^{KIS2^2Ki}'' 

[ '^{KlS2^2Ki} R2 
Gk < G(k-^^ iteratively for k = 3, . . . ,n — 1 and 



S2S1, K^2) = [ S2K1 K2 ] and R 



Repeating the computation for G 
writing X(fc) 



(fe) 



,xD^and 



r(fc-l) 



(k)X[k), 



fc, we arrive at the desired result with G = G(^n), S = S(n), K = K(^n) 



at each iteration 
and R = R(^n) as 
stated in the lemma. 

To see that R + Q{K''K} is lower 2x2 block triangular, we note that K'^K may be 
expressed as follows: 



K^K 



kIKi 
KIS2K1 



k\sIk2 
kIk2 



Klsl^2K2 



KnSn<^2Kl KnS„ 



iK2 



n-1 



^l'^n^2^n 



Note that since K'^K is by definition a Hermitian matrix, the 2x2 block elements above 
the diagonal blocks are the Hermitian transpose of the corresponding elements below the 
diagonal blocks. It follows therefore that the imaginary part of the block {K^K)jk at 
block row j and block column k must satisfy the relation: ^{[K'^ K) jk} = —'^{{K''K)kj}'^- 
However, from the expression for R derived above and its symmetry, we already have that 
if A; > j: 

R,k = Rl = ^{KlSk^.^iK,}^ = Q{iK^K)k,V. 

Therefore, the off-diagonal upper block elements of R cancel those of Q{K'^K} when they 
are summed and we conclude that the matrix R + Q{K'' K} is a lower 2x2 block triangular 
matrix. 

Recall again the partitioning of R as R = [Rjk]j,k=i,...,n with Rjk G M^^^ and of K as 
K = [ Ki K2 ... Kn ] with Kk E C"^^^. We may now state the following result: 

Theorem 4 A linear quantum stochastic system G = {S, Kx, ^x^Rx) with n degrees of 
freedom is realizable by a pure cascade of n one degree of freedom harmonic oscillators 
(without a direct interaction Hamiltonian) if and only if the A matrix given by A = 2Q{R + 
'is{K^K}) is a lower block triangular matrix with blocks of size 2x2. // this condition is 
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As we already saw in the proof of Lemma 



whenever k > j and Rkj = '^{K^.Kk}'^ 



satisfied then G can be explicitly constructed as the cascade connection Gn<Gn~i< ■ . .<Gi 
with Gi = {S, KiXi, ^xj RiiXi) , and G^ = (/, K^Xk, ^x^Rkk^k) for k = 2, . . . ,n. 

Proof. The proof of the only if part follows directly from Lemma [3} as follows. If G can 
be realized by a pure cascade connection of n one degree of freedom harmonic oscillators 
then by the lemma, R + Q{K'^K} is a lower 2x2 block triangular matrix. However, since 
is 2 X 2 block diagonal, it follows that the matrix A = 2Q[R + Q{K'^K}) is also a lower 
2x2 block triangular matrix. 

Conversely, the if part of the proof can be shown by explicitly constructing a pure 
cascade connection of n one degree oscillators that realizes G. If the A matrix associated 
with G is lower 2x2 block triangular then so is the matrix ^Q~^A = —^QA = R+Q{K'^K}. 

this structure implies that Rji^ = Q{KjKk} 
k < j. Now, using the notation of Lemma 
[3| let us define the one degree of freedom harmonic oscillators G^ for k = l,...,n as 
Gi = {S, KiXi, ^xjRiiXi), and Gk = (/, KkXk, ^x]^RkkXk) for k = 2, . . . ,n. It follows from 
Lemma |3] that Gn < Gn-i < ■ ■ ■ < Gi = {S, Kx, ^x^Rx). That is, this cascade connection 
realizes G. 

Theorem|4]has a direct consequence on the weaker notion of transfer function realization 
of a linear quantum system. The main result is the following corollary: 

Corollary 5 A linear quantum system G = {S,Kx, \x^Rx) is transfer function realizable 
by a pure cascade connection of one degree of freedom harmonic oscillators if and only 
if there is a symplectic transformation matrix V such that the linear quantum stochastic 
system G' = {S, KV~^x, \x'^V~^ RV~^x) has an A matrix which is lower 2x2 block 
triangular. 

Proof. By Definition [T| G is transfer function realizable by a pure cascade connection if 
and only if there exists a symplectic matrix V such that G' = {S, KV~^x, ^x^'^V"^ RV~^x) 
is realizable by a pure cascade connection. But from Theorem |4] this is true if and only 
if the A matrix associated with G' (i.e., A = 2Q{R' + '=s{K'''K'}) is lower 2x2 block 
triangular. 



5 Passive linear quantum stochastic systems 

In this section it will be shown that the class of passive linear quantum stochastic systems 
(as defined below) are transfer function realizable by a cascade connection. In [TU] it 
has been shown by a constructive algorithm that a "generic" sub-class of such systems are 
transfer function realizable by pure cascading, the generic systems being required to satisfy 
assumptions on the distinctness of the eigenvalues and invertibility of certain matrices. In 
this section we remove such assumptions, and show by exploiting the algebraic structure 
of passive systems that the result is valid in general for all passive systems. 
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For A; = 1, . . . , n, let flfc = {qu + ipk)/'^ be the annihilation operators for mode 



fine a = (ai, . . . , a„)"^. Then a satisfies the CCR 
diag(/„, — /„). Moreover, note that {a^,a^)'^ = [ E"'' with 



a 



k and de- 

T 

at 



2 2^ 


... 




1 

2 

















We also make note that 
/ we have the identities 



E 

E# 



2 [ E^ E"^ ] and from the relation 



E 

E# 



2 [ Et E^ 



EEt = J/2 = E^E"^; EE'^ = = E#Et. 



(2) 



Therefore, we also have 





E 


-1 


a 


X = 


. ^* . 







2 [ Et E^ 



A system G = {S, Kx, ^x^Rx) is said to be passive if we can write H = \x^ Rx = ^( 
and L = Kx = Ka for some complex n x n Hermitian matrix R , a complex m x n 
(here m again denotes the number of input and output fields in and out of G) matrix K, 
and some real constant c. As discussed in here the term passive for such systems is 
physically motivated since they can be implemented using only passive components like 
optical cavities, mirrors, beamsplitters and phase shifters; this follows from Theorem 5.1 
of [7] and the constructions shown in section 6 of that paper. Also shown in fO], we can 
express ^a^Ra and Ka in the form ^a^Ra = ^X"^3ft{Eti?E}a; — | J2]=i ^jj — KT,x. 

Therefore, we may set R = 3f?{Et^E} and K = KH. Note also from [9] that the 2 x 2 



block diagonal elements {Rjf j 



n} is of the form Rjj = Aj/2 for some Xj e M for 



all j. Now we shall derive some properties of A and show that there exists a unitary and 
symplectic matrix that transforms it into a lower 2x2 block triangular matrix. 

Lemma 6 [ E^ Et fA [ Et E^ ] = diag{M,M*), where M = iEeEt(^ - ii^t^). 

Proof. For the proof, we exploit the identities ^ as well as the following easily verified 
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identities: EGS-^ = = S*0S"''. Using these identities, we have the following: 

= 2[ St ]^e(3R{Sti?S} + s^isti^-ti^s}) [ St s^ ] 
= 2[s^ St ]^eQ(sti?s + s^i?#s#) 

-^(Sti^tKS - S^if^ir#S#)^ [St s^ ] 

= 2[s^ St ]^e[ isti?-istirti^ Ij:^r# + {j:'^k^k*] 
= diagQsest^- ^sesti^ti^, 

-s#es^i?# + ^s#es^i^'^i^# 

Then we have the following theorem: 

Theorem 7 Let U be the complex unitary matrix in a Schur decomposition of the matrix 
M of Lemma @- M = WMU, where M is a lower triangular matrix. Then the matrix 

y = 2[ St S^ ]diag{U,U*)[ S^ St f 

is a real, unitary, and symplectic matrix that transforms A into a lower 2x2 block trian- 
gular matrix: VAV'^ = A, where A is a real lower 2x2 block triangular matrix. There- 
fore, every passive linear quantum system has a transfer function realization by pure cas- 
cading and such a realization is obtained by applying the construction of Theorem to 
G' = {S,KV^x,^x^VRV^x). Moreover, each of the one degree of freedom oscillator in 
the cascade will also be passive. 

Proof. The existence of U is guaranteed by the well known result that every complex 
matrix M has a Schur decomposition of the form M = WMU with M lower triangular. 
Note then that we also have M* = U'^M'^U'^ . Let V be as defined in the theorem. Then 
by Lemma |6] the following is true: 

VAV^ = 4[i:^ S'^ ]diag(M,M#)[ S'^ St = 4(StMS + S^M#S#) = 83f?{StMS}. 

Now, since M is a lower triangular matrix, it follows by inspection (using the special struc- 
ture of S) that A = 83fJ{StMS} is a lower 2x2 block diagonal matrix, as claimed. That 
V is real follows from the fact that we may write V = 2(Stf/S + S^f/#S#) = 43fJ{Stt/S}. 

That it is unitary follows from the observation that \/2[ St S"^ ] and \/2 



S 

s# 



are uni- 



tary (as a consequence of (|2|) and that diag(f/, f/*) is also unitary. To see that V is also 
symplectic define b = Ua and z = Vx. By the unitarity of U we have that b and 6* satisfy 
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the same the commutation relations as a and a* (i.e., h is again an annihilation operator). 
Then we have 



Vx = V2[ St 



2[ St S^ ]diag(f/, f/#) 



[ St S^ 



b 

b* 



X] 



Now, this implies that z consists of the canonical position and momentum operators as- 
sociated with the modes in b and satisfies the same CCR as x. But since z = Vx and V 
is real, preservation of the CCR implies that V is necessarily a symplectic matrix (this is 
standard knowledge in quantum mechanics; see, e.g., [2H1 section III]). 

Using the fact that V^^ = V'^ = V^t established above, it follows from Theorem |4] that 
the passive quantum system G is transfer function equivalent to G' = {S, KV^x, ^x^VRV^ 
whose A matrix is lower 2x2 block triangular. Let K' = KV^ = [ K[ K'2 . . . K'^ ] 
and E! = VRV'^ = [R'jk\j,k=i,...,n- By Theorem |4 we have G' = Gn< G^-i <...<Gi with 
Gk = {Sk, K'l^Xk, \x]:R'kk^k), Si = S and Sk = I for k > 1. We now show that each Gk is 
passive. Recall that K = and write K = [ Ki ... ] with Kk G C™^^ Using ^ 
we have that K'x = KV'^x = KWa. By expanding both sides of the equality K'x = KWa 
and collecting and equating terms of the same index, it follows that K'g^Xk = {KU'^)ka,k, 
where {KU'^)k is the /c-th C"*^^ block component of KU'^ . On the other hand, since G is 
passive we have that Rkk = for some \k e M, and recalling that R = 3fJ{Sti?S}, it 



follows by inspection (after some algebraic manipulations using ([2])) that i?'^^ = X'f^l2 for 
some A';. G M. Thus, x'j^R'f^^Xk = X'l^alak + c for some constant c and we conclude that each 
Gk is also passive. 



Example 8 Let G = (/, Ka, ^a'^Ra + |) be a passive system with R 



K 
R 



. Here K 



1 + 0.51 -2 + i 
-5 - 2i 3 - 4i 

0.5/2 0.25(/2 - J) 
0.25(l2 + J) 0.75/2 

-0.6933 + 0.0039i 0.2244 - 0.68491 
0.7204 + 0.02091 0.2312 -0.65361 

Then by the formula of Theorem we have 



0.5 + 0.251 
-2.5-1 



-0.25 + 0.51 
1 - 2.51 



2 1 + 1 
1-1 3 
-1 + 0.51 -0.5-1' 
1.5-21 2 + 1.51 



By Lemma 
and M = 



6 and Theorem 



and 
and 



-14.8390 -0.79121 
0.6344 - 0.22251 



?, we have that U 





-0.2235 - 0.45881 



V 



-0.6933 0.0039 0.2244 0.6849 

-0.0039 -0.6933 -0.6849 0.2244 

0.7204 -0.0209 0.2312 0.6536 

0.0209 0.7204 -0.6536 0.2312 



Let K' = KV-^ 
K' = 



[K[ K'^], then 



-0.9144 - 0.74411 0.7441 - 0.91441 
3.4433 + 1.26211 -1.2621 + 3.44331 



-0.1926 -0.36841 0.3684-0.19261 
-0.1679 -0.15011 0.1501 - 0.16791 
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and 



R' = VRV-^ 
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12) 



Rl2 
-^22 



0.7912 

0.7912 

0.1113 -0.3172 

0.3172 0.1113 



0.1113 0.3172 

-0.3172 0.1113 
0.4588 

0.4588 



Therefore, by the theorem, G is transfer function realizable by G' = (/, K'x 
easily checked that R' + '^{K'K} is lower 2x2 block triangular: 



It is 



R' + "^{K'^K'} 



0.7912 14.8390 

-14.8390 0.7912 

0.2225 -0.6344 0.4588 0.2235 

0.6344 -0.2225 -0.2235 0.4588 



therefore G' can be realized by a pure cascade connection of one degree of freedom harmonic 
oscillators. According to Theorem^ G' = G2 < Gi with Gi = {I,K[xi,^xjR'iiXi) and 
G2 = {I , K2X2,\x'^R'22X2) ■ It is casHy inspected that both Gi and G2 are passive. A 
quantum optical realization of G' is illustrated in Fig. [1} 




Figure 1: Realization of G' as the cascade connection of Gi and G2- Gi and G2 are each 
reahzed by an optical cavity and a phase shifter; see, e.g., [221 1231 [Z| for a discussion of 
these devices. Dark rectangles depict fully reflecting mirrors, while light rectangles depict 
partially transmitting mirrors; an optical cavity is formed by bouncing light back and 
forth between two mirrors. Here 9n = —2.4585, 6^12 = 0.3513, 6*21 = —2.0525, and 6*22 = 
—2.4121, and the partially transmitting mirror Mjk in the optical cavity Gj, j,k = 1,2, 
have the coupling coefficients 711 = 1.3898, 712 = 13.4492, 721 = 0.1728 and 722 = 0.0507, 
respectively. The resonance frequencies of the optical cavities of Gi and G2 have a detuning 
of 0.7912 — 0.4588 = 0.3324, with Gi having the higher resonance frequency. 
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6 Conclusions 



We have derived a characterization of hnear quantum stochastic systems that can be real- 
ized, in a strict or transfer function sense, by a cascade connection of one degree of freedom 
quantum oscillators alone, without requiring any direct bilinear interaction Hamiltonian 
between these oscillators. The results are constructive in that if a system can be realized 
by a cascade connection, it is explicitly shown how to construct it. Then it was shown 
that the sub-class of passive linear quantum stochastic systems is always transfer function 
realizable by a pure cascade connection. 
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